arXiv: 1503.00127v2 [physics.soc-ph] 1 Apr 2016 


i 


How crude oil prices shape the 
global division of labour 

Francesco Picciolo 1 ’*, Andreas Papandreou 2 , Klaus 
Hubacek 3 & Franco Ruzzenenti 4,5 
1 Department of Chemistry, University of Siena, Via Aldo 
Moro 2, Siena, 53100, Italy 
department of Economy, University of Athens, 8 
Pesmazoglou Street, Athens, 10559, Greece 
department of Geographical Sciences, University of 
Maryland, CollegePark, Maryland, USA 
department of Economics and Statistics, University of 
Siena, Siena, Italy 

5 Institute of Sociology, Jagiellonian University, Krakow, 
Poland 

*To whom correspondence should be addressed; E-mail: 
picciolo.fr@gmail.com 

INTRODUCTION 

In the aftermath of the oil crisis of the early 1970s, the re¬ 
lationship between oil prices and economic growth became a 
focal point of the scientific discourse and public debate. In 
1983, James Hamilton published an influential article show¬ 
ing that an oil price increase had preceded all but one reces¬ 
sion in the United States since the end of World Il ETIl . Since 
then, a large number of empirical studies have looked into the 
connection between oil prices and real economic growth and 
frequently found a significant negative correlation lfTll27l . The 
importance of this link between oil price and economic growth 
was less clear after the second oil shock[[21 E3] ESI E23 - Recent 
studies, with more refined statistical tools and price specifica¬ 
tions, have accomplished in restoring the between the oil price 
and economic growth lfl4l l22l 1331 l35l 36). There is now a gen¬ 
eral consensus on the notion that this did not cease but it has 
become more complex in terms of direction (anticyclical and 
procyclical), typology of shocks (demand or supply) and lag 
patterns lf9| l28l l33l [351 . This line of research tried to explain 
this tight relationship, given that the cost of energy is but only 
a small part of GDPf9) but satisfactory explanations have re¬ 
mained elusive j9] [27 '28, 371. Interestingly, this research on 
the oil price and economy economic indicators seems to have 
entirely ignored the transport sector, which is heavily reliant 
on refined crude oil products, and its role in shaping the global 
division of labour. In the post-war period, world trade grew at 
a faster pace than world GDP ||25) . According to recent studies 
on globalization, the remarkably high rate was propelled by a 
dramatic decline in international transport costs l5l IT3l [131 l25l , 
Perhaps, the notion that trade grew amid globalization because 
of transport should not come as a surprise. What is more 
surprising, but is closely related or even a corollary, is the 
fact, that intermediate and capital goods, in the last decades, 
grew faster than final products and now account for the largest 
part of trade in OECD countries ED- While, in the after- 
math of World War II international trade mainly concerned fi¬ 
nal products, the second wave of globalization (since the late 
1980s) extended to intermediate products and capital goods. 


and the integration of factors market as another important 
effect j7). This process led to the fragmentation of production 
internationally |26l . Disregarding the transport sector, most of 
the scholars focused their attention on other factors in order 
to explain the fragmentation of the global value chain, like 
the pursuit of cheap labour or more favourable environmental 
legislations If7l [26] [29 1. Amador and Cabral recently suggested 
that the strong increase of trade associated with the develop¬ 
ment of the global value chains (GVCs) in the 1990s coincides 
with a period of low oil prices, although admitting that there 
is little empirical evidence linking these two factors (3). These 
findings suggests the importance of assessing the impact of oil 
shocks in an internationally integrated system rather than on 
a national base only. Furthermore, the new issues posed by 
climate change demand a deeper understanding of the nexus 
between energy consumption and the global economic struc¬ 
ture. Our study addresses the connection between oil price 
and the global economy, by means of network theory and 
Markov chain theory, with the aim of understanding how the 
GVCs expanded and shrank following price changes in crude 
oil, between 1960 and 2011. In contrast to previous analyses, 
that progressed by refining price specifications and statistical 
methods, we observed the correlation of economy with crude 
oil price (Brent), but we changed the macro-economic vari¬ 
ables under investigation. We first applied network theory to 
trade imbalance and bilateral trade to understand how these 
two global measures of trade are linked to the oil price. These 
two quantities are thereby used to introduce the cycling index 
that builds on Markov chain analysis to assess the amount of 
value that is conserved across direct and indirect relationships 
in trade. With this measure, we looked at the share of cyclical 
value the share of value that returns to the starting point, along 
different paths in the world trade network. 


ANALYSIS AND RESULTS 
A. Balance of trade per country: trade (im)balance 

The balance of trade is the difference in value of exports and 
imports (see methods). Many have viewed the existence of 
large current account imbalances between large economies as 
a possible cause of the financial crisis |[T2ll40l . There is grow¬ 
ing evidence that current account (im)balances are correlated 
to oil prices worldwide 140) . The reason for this correlation 
lies in the burden placed on imports (or exports, for oil ex¬ 
porting countries) by energy commodities, but also in mone¬ 
tary policies aimed at regulating inflation (which is correlated 
to oil price) P0) . The analysis has been performed on a yearly 
basis between 1960 and 2011, on aggregate trade flows (total 
import/export for every country) of all the reporting countries 
in the world. Data are taken from Gleditschs lfl9) and BACI 
datasets lfT8l . The fluctuation in the balance of trade and the 
variation of the adjusted crude oil price are slightly, yet sta¬ 
tistically significantly, negatively correlated: the linear corre¬ 
lation coefficient is -0.32 (see tab.l). It is noteworthy that in 
a network where flows tend to be balanced at every vertex, 
the matrix tends to be symmetrical (which means that entries 
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in the upper triangular matrix mirror those in the lower tri¬ 
angular one) l46l . In other words, symmetric weights (flows) 
between every pair of vertices is statistically the simplest way 
to balance ingoing and outgoing flows at every vertex. A local 
symmetry (exports equals imports) tends to produce a global 
symmetry (export from i to j equal export from j to i). We 
thus expect that the balance of bilateral trade in the world trade 
web (WTW) to be negatively correlated to oil prices because 
we observe a negative correlation of the balance of trade lo¬ 
cally. It should be noted though, that this is just statistical 
relationship, obtained by imposing the local balance as a con¬ 
straint in the null model fl46ll . 


B. Balance of bilateral trade: trade reciprocity 

A global measure for evaluating the balance of trade be¬ 
tween every pair of countries is the weighted rccinrocity |[46i . 
Reciprocity is a first-order property, meaning that it concerns 
the direct relationship of nodes with the nearest topological 
neighbours (one link-length). Reciprocity has proven to be 
an helpful measure in understanding the effects of the struc¬ 
ture on dynamic processes, explaining patterns of growth in 
out-of-equilibrium networks lfl7l l52l . and starting to evaluate 
higher order properties |30jS307 3131511 • The reciprocity for 
weighted networks, r w , is defined as follows: 

w EiEj min Kj.«'3i] 
r Ei £,■ wn ’ 

If all flows are perfectly reciprocated/balanced then r w = 
1. If they are unidirectional r w = 0. The correlation of 
r w with oil price is -0.70. As the oil price increases (de¬ 
creases), reciprocity decreases (increases). This result shows 
that changes in the international bilateral trade stmcture and 
the oil price are more intimately and inversely linked than 
expected by observing only global volumes or distances of 
trade ll5l fl3l i25l . These results deliver two important pieces of 
information:!) given that the correlation of the oil price with 
reciprocity is higher than the correlation with imbalance, the 
nexus between oil price and reciprocity (bilateral trade bal¬ 
ance) cannot be reduced to the correlation between oil price 
and imbalance, despite the expected symmetrical effect (see 
above); 2) this tight correlation is not explainable with trading 
agreements, as it is ubiquitous l42l . nor with a general abate¬ 
ment of barriers Il25l . as it has a discontinuous trend in time 
(see SI). Can we link reciprocity to some structural effect in 
production, like the development of GVCs? As previously 
stated, reciprocity is a first-order property, though, unless we 
assume that production chains involve only two sites in a row, 
GVCs should be pertinent to higher-order properties of net¬ 
works. In the SI we show (fig. 4), with a heuristic model 
based on three countries, one product and two factors, how 
reciprocity increases when the production chain expands, in¬ 
volving second order properties of the network (neighbours 
of neighbours, or indirect relationship). Within this simple 
model, the shift from a single-country production chain to a 


multi-country production chain will always increase the reci¬ 
procity of the network, independently on the distribution and 
share of the total volume traded. We thus expect that shift¬ 
ing the production sites abroad increases the reciprocity of the 
network. However, to check this hypothesis we must extend 
the analysis to longer paths of the productive chain, involving 
more than three nodes (like in the heuristic model) and, most 
importantly, encompassing indirect relationships. 


C. Markov chain analysis: cyclic paths of value in the world 
economy. 


The largest share of trade in the world economy involves 
inputs to production (raw materials, intermediate and cap¬ 
ital goods) f3Tj. In the modern economy countries import 
these production inputs and export final products or inter¬ 
mediate goods that are further processed elsewhere often in¬ 
volving numerous production stages in many different coun¬ 
tries. At every step in the global production chain value added 
is embodiedf3|. We are interested in detecting the share of 
traded value that is conserved throughout the stages of GVCs, 
or the initial value at the beginning of a global production 
chain, like mass particles in ecological networks that are con¬ 
served throughout every stage of a food chain fl6l (see meth¬ 
ods for a detailed description of the concept). Thus, we want 
to assess the share of trade that is cyclical that runs a cyclic 
path, along paths of a given length S, in the WTW. By means 
of Markov chain theory, we can statistically evaluate the prob¬ 
ability of an ’’elementary” trade, i.e. the amount of value em¬ 
bodied in raw materials or intermediates that is conserved in 
a product, going from country i to country j and return to the 
initial country throughout all the possible direct and indirect 
trade relationships lITOl [34l . For example, Fig.l illustrates all 
the possible cyclic paths that a particle (a unit value of trade, 
in our case) can follow from country i within the first 4 steps 
(i.e. stages of production). 

In order to assess the share of cyclical value over the total 
value traded, we need a normalized measure, like the previ¬ 
ous measures of imbalance and reciprocity. We indicate with 
the cycling index F' 5 - the share of trade that comes back to 
the starting country in S steps with S = 2,, oo (see methods). 
In Fig. 2 trends of the cyclical quantity are shown for the 
WTW, for S = 2,3,4, oo. We observe that the percentage of 
cyclical value inside the network and oil price are negatively 
correlated, at various degrees. The lower the number of steps 
(i.e. the shorter the production chain) taken into account while 
evaluating the higher the negative correlation. The corre¬ 
lation from F^ to r(°°) goes from -0.85 to -0.62, exhibiting 
a much higher score compared to imbalance and reciprocity 
(see tab.l). 
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FIG. 1. Cyclic paths up to four steps. The possible cyclic paths that 
a trade can follow on a network within 4 steps. Node i (represented 
by the black dot) is the starting and ending country of each path. 
T (S) evaluates the share of the total trade that follow a cyclic path of 
length up to S. 


DISCUSSION 

D. Digging into the tight relationship between cycling and oil 

price. 

One may argue that a general increase in exports and im¬ 
ports driven by global GDP growth, would inevitably lead to 
raise the share of mutual trade. In this view, the correlation 
between cycling and oil price could be brought back to the 
GDP-oil price nexus. Moreover, it is also possible that oil 
prices inflating commodity prices and automatically increas¬ 
ing nominal trade flows, could explain this tight correlation. 
In this latter case, the relationship would boil down to the link 
between inflation and oil price. Thirdly, it is possible that the 
correlation between oil price and cycling is a multiplicative 
effect (re-spending effect of petrodollars) of the correlation 
between imbalance and oil price, linking oil producing with 
oil consuming countries. This is conceivable on the notion 
that oil exporters tend to import goods and services from oil 
importers. In summary, all three may be included in the no¬ 
tion of first order properties of the network contrary to cycling 
that is a second order property. We can show, using a Null 
Model as a benchmark that the negative correlation between 
the cyclical value inside of the network and oil price does not 
depend on first order properties (direct relationships), but de¬ 
pends on higher-order properties of the network (indirect re¬ 
lationships). The chosen null model is the Weighted Directed 
Configuration Model (WDCM). The WDCM is a well-known 
and utilized null model that preserves the quantity of trade 
and the distribution of trade relationship of each country ll46l . 
We performed the analysis of the null model and we evalu¬ 
ated the correlation of the expected cyclical value with the oil 
price (see methods). The correlation between the expected 
T^’s and the oil price are significantly lower than the ob¬ 
served ones. They range from -0.32 to -0.38. Notably by 



Year 

FIG. 2. The trend of the cyclical flow index, T^ s \ of the World Trade 
Web calculated in eq.(5) from 1960 to 2011. From lighter to darker 
color: T (2) (orange line), T® (red line), T^ (purple line), and T (o °^ 
(black line). 


using the WDCM we preserve: 1) the nominal value of all the 
exports/imports sequence; 2) the global value of trade; 3) the 
difference between import and export (balance of payments) 
at every vertex. The comparison of expected correlation from 
the Null Model with the observed correlation assures us that 
the latter is not due to the global volume traded, nor to the 
nominal values of trading relationships or the distribution of 
trade imbalance: it depends on the specific architecture of 
trade flows globally, which can only be understood by the way 
GVCs unfold. In order to understand the underlying process, 
it is also instructive to observe the trend of cycling indexes 
of different path-length over time. In Fig.2 and Fig.3, the 2- 
steps length peak of the early 1990s is lower than the one of 
the late 1960s, meaning that the higher orders contribution to 
the cycling has become more and more important, following 
the decreasing trend of the 1970s. If we look at the trend of 
the normalized values of the two, three and four steps cycling 
compared to oil price (Fig.3), it is evident that until the mid- 
1970s the three degree of cycling overlap, while between the 
first and the last oil crisis (2008) the curves of the three and 
four step cycling stand above that of the two step cycling. The 
longer-than-two step cycling also displays higher growth rates 
until the peak of the 1990s, suggesting that the ongoing glob¬ 
alization was characterized by longer paths of cycling. Longer 
cycles in the value chain probably underline the growing share 
of intermediate goods in trade and a process that led to a more 
interdependent global economy (see SI). 


E. The second wave of globalization and the role of transport 

We hypothesize that the transmission mechanism behind 
the correlation between oil price and GVCs, statistically mea¬ 
sured by cycling, is the transport sector. Oil prices impact 
on transport costs, making international outsourcing more or 
less profitable fiTi l48l . If this hypothesis is correct, we ex¬ 
pect oil prices to impact on cycling, by influencing the length 
of GVCs. We tested Granger causality, between oil price and 










4 


the three network measures here considered (imbalance, reci¬ 
procity and 2-step-cycling), in both directions, with one lag 
specification and a significance level of 5% fl20l . Only in the 
case of cycling we can reject the null hypothesis that oil price 
does not Granger cause 2-step-cycling. The Granger test, 
however, indicates that causation runs in both directions. In¬ 
deed, the most important information brought about Granger 
causality test is that 2-step-cycling is cointegrated with oil 
price, meaning that these two variables follow the same trend, 
pointing to a long-term dynamics. By enlarging the scope 
of the analysis from first order properties of the network to 
higher order properties, we obtained not only a higher corre¬ 
lation to oil, but also a cointegration, meaning that the analysis 
through the cycling quantity are able to detect and capture the 
global changes of the pattern of WTW. It thus seems plausible 
to believe that transport is the nexus between GVCs and oil 
price, as we expect changes in the structure of production to 
occur in the long-term. However, there are many other factors 
that influence the global division of labour that would need 
to be included to better explain this relationship. What seems 
somewhat surprising, given the declining costs of transport, 
is the fact that the process of international integration peaked 
in the 1990s and declined in the 2000s. Indeed, the second 
half of the 1990s witnessed the onset of the second wave of 
globalizationQ. This apparent contradiction is probably ex¬ 
plained by observing single-country cycling (the portion of 
cycling passing through a single vertex, see fig. 4). Dur¬ 
ing the 2000s the single-cycling of fast-growing economies, 
like China, increased dramatically, climbing the ranking of 
the World’s economies. While the cycling index of developed 
economies remained constant, China showed a rapid growth 
in the share of cyclical value of its trade, meaning that a large 
part of the GVCs began passing through this country. Some 
economies like China attracted a significant portion of the cy¬ 
cling value across the World, concentrating the flows into a 
few hubs. The decrease in cycling globally is consistent with 
the emergence of hubs, which are topologically like stars (see 
fig.4 in SI). It is plausible that outsourcing initially, from the 
late 1980s, was propelled by road transport, over medium dis¬ 
tances and involving many countries (inflating global cycling) 
whereas in the second stage, from the 1990s, was fuelled by 
air and cargo shipping, concerning longer distances and few 
hubs (reducing global cycling). According to IEA, interna¬ 
tional cargo shipping, mostly in non-OECD countries, dis¬ 
played the highest growth rates among the different transport 
modes between 1990-2000 1481 . It is noteworthy that the first 
major shift coincides with a sudden leap in efficiency in road 
transports triggered by the oil crisis iBTll . whereas the second 
change followed a dramatic decrease in the international trans¬ 
ports costs, both in air and cargo shipping, following the in¬ 
troduction of a more efficient aircraft fleet and the containers 
svstem |[25l . To investigate this hypothesis, we assessed how 
the role of distances in shaping the cyclical value has changed 
over time. If we trim the network at different distance thresh¬ 
olds (removing all the links placed between a couple of ver¬ 
tices above a certain distance value) and calculate the cycling, 
we observe that there have been a bifurcation between short 
and long distance cycling in the early 1990s (see fig.6 in SI). 



Year 

FIG. 3. The normalized crude oil price (black solid line with sub¬ 
tended grey area) is shown. The trend of the normalized value of 
T (2) (black solid line), T(3) (dashed line), T^ (dotted line) are also 
shown. The linear normalization is done based on the maximum and 
minimum value reached by each quantity. The final scaled values lie 
between 0 and 1. 


While before the 1990, the cycling indexes trimmed at dif¬ 
ferent distances were moving coherently, after the 1990s the 
share of short-distance cycling index (less than 2500 km) be¬ 
gan to decline compared the share of long-distance cycling in¬ 
dex, which kept growing (more than 2500 km). This seems to 
indicate that the GVC’s shifted from a regional domain (clus¬ 
ters of neighbouring countries) to a global (oversea relation¬ 
ships) domain. 



R 

p-value 

b t 

-0.32 

0.02 

r w 

-0.70 

< 10 3 

p( 2 ) 

-0.85 

< 10 3 

p(3) 

-0.83 

< 10 3 

p( 4 ) 

-0.82 

< 10 3 

p(oo) 

-0.62 

< 10 3 


TABLE I. The linear correlation index between b t , r w , T (S ^ (for 
S £ {2, 3,4, <x)}) and the crude oil price are reported, in the third 
column the respective p-values are shown (95% confidence). 


CONCLUDING REMARKS 

Global trade, in the age of the second globalization, 
has entangled national economies, by interconnecting pro¬ 
duction sites internationally, in a fashion that is still 
underappreciated(3] |7) [20,23]. This historical process of 
vertical integration has been producing global value chains 
(GVCs) wherein goods travelling across countries augment 
their value at every stage of production^. We investigated 
this cyclic path of value across countries by means of network 
theory and Markov chain theory iflOl [34] [46] |. In this way, we 
show that over a longer time period the oil price has a striking 
correlation with the structure of trade globally. A correlation 
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(b) 



(c) <d) 


FIG. 4. The ranking position for seven countries according to the 
value of (panel a) and 7 ; - oc ' 1 (panel b), T'l' 2 ' 1 (panel c) and r^°°^ 
(panel d). The data range covers from 1960 to 2011. The seven 
countries are the following: China (orange), Germany (red), France 
(bleu), Great Britain (purple), Japan (green), Russia (cyan), USA 
(black). We chose these countries for their interesting patterns in 
order to reveal how the variables 7 { S> and T^ can explain struc¬ 
tural changes in economy. Russia position is always higher in the 
relative than in the absolute rank, and for S = 2 than for S' = oo. 
Moreover we can observe a decline in both starting in late 80s and a 
rise since around 2000. Results seem in line with the USSR disso¬ 
lution, explained by a Russian chain value made up substantially by 
bilateral exchanges with other former Soviet Republics. A remark 
on the performance of China and Germany: the former has been on 
the rise in the last decades and we point out the pervasiveness of its 
chain value in absolute value reaching the top of the rank for in 
the last years of our analysis. In the same period, Germany overtakes 
USA share of cyclic trade (S = oo ), as emerges from panel d. 


that increases with the scope of the analysis, from first or¬ 
der properties of network (one link distance), to higher order 
properties. The worldwide sum of country trade imbalances 
show a weak correlation of -0.32. The correlation increases by 
engulfing bilateral relationships between countries on a global 
scale (reciprocity, -0.70). Finally, the highest correlation, up 
to -0.85, is observed when we involve more complex pattern 
at the global level. By means of statistical mechanics of net¬ 
works (exponential random graphs) we were able to demon¬ 
strate that this remarkable correlation can only be explained 
with higher-than-one order properties of the network, indicat¬ 
ing that GVCs and structure of trade are intimately linked to 
oil price. We hypothesize that this tight relationship points 
to the role of transports in determining, in the long run, the 
extent and the way production sites connect internationally. 
By looking closely at the single-country cycling index (fig. 
4) and by dissecting the global cycling according to different 
distance thresholds (fig. 1 in SI), we individuated two struc¬ 
tural breaks and two phases of the second wave of globaliza¬ 
tion (the second unbundling). The first, which started in the 
1980s and peak in the 1990s, was featured by shorter distance 
and increasing cycling, the second, that probably lasted until 
the economic crisis, was featured by longer distances, declin¬ 


ing global cycling and increasing cycling of China (star-like 
structure). Our results suggest that the transmission mecha¬ 
nism between oil price and economic growth lays not only in 
the labour or retail markets (via inflation), but also, more pro¬ 
foundly, in the structure of production globally. Furthermore, 
in a more general perspective, our results indicate that the pro¬ 
duction structure could be approached as an energy system, 
constrained by energy efficiency in the transport sector. This 
view of the economic system builds on the work of scholars 
like Ayres (6), for approaching growth as a product of energy 
efficiency, but also on fundamental advances in the study of 
allometric scaling, aimed at explaining the structure and size 
of many biological processes as the result of general features 
of efficient transportation networks®. Nevertheless, in order 
to establish a clear causation among factors, further research 
should tackle two aspects more in depth: 1) the dynamic pro¬ 
cess between input/output matrices of national economies and 
sectoral trade; 2) the evolution of energy efficiency and costs 
of transports globally, in the long run. The authors of this ar¬ 
ticle believe this is crucial for understanding the role of oil in 
the present economic system, given that this is not yet a fun¬ 
gible source of energy in the transport sector, and for paving 
the way for a prosperous economy freed from fossil fuels. 


Appendix A: Materials and methods 
1. Description of the dataset. 

In the following sections, a brief description of the analysed 
networks is given. We analyse the series of yearly bilateral 
data on exports and imports among world countries from the 
Gleditsch’s database fT91 . from 1960 to 1997. From 1997 to 
2011, we employed international trade data provided by the 
BACI database lfl8l . All data are in millions of current U.S. 
dollars and all they are freely available. Also the crude oil 
price data are freely available (TT). 


2. Evaluating the balance of trade. 

In network approach the World trade web is described by 
matrix IF, where each node i is a country and w l3 is the trade 
from country i to country j, we can evaluate trade imbalance 
as follows: 


6 * = 


E a min[ S r,sn 
E* Ej w ij 


(Al) 


where s) n = ]TT Wji is the in-strength (total import) of 
node i and s° ut = ^ ■ is the out-strength (total export) 
of node i. In doing so, we do not distinguish between trade 
surplus or deficit (positive or negative value), instead we cal¬ 
culate only the amount of trade that is balanced between im¬ 
ports and exports for every country (i.e. the minimum of total 
imports and exports); 6* = 1 if exports and import are equal 
for all the countries, 0 > b t > 1 otherwise. 
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3. Evaluating patterns of traded value. 


We want to assess the probability of a percentage of trade to 
be in a country j, starting from i, after a number of production 
steps. Markov chain analysis allows us to statistically deter¬ 
mine this. The probability is calculated according to the trade 
relationships between country in WTW reported bvll8l fl9l . 

Calculating transition probabilities. In probability theory 
a Markov chain is a stochastic process defined on a discrete 
state space satisfying the Markov property. It is a set of ran¬ 
dom variables, representing the evolution of a certain system, 
without memory: each actual state of the system just depends 
on the previous one. The changes of states are called transi¬ 
tions and the Markov chain can be described by a Markov ma¬ 
trix, M, whose elements m VJ represent the probabilities from 
a state i to a state j (transition probabilities or transition rules). 
Given the memory-less nature of a Markov chain it is not pos¬ 
sible to predict the state of a system, in a given time t, but 
it is possible to predict its statistical properties. In this paper 
the transition probabilities are set according to the elements 
of the matrix W that describes the World Trade Web, i.e. the 
trade network. A link, w. rj , is assigned to any import/export 
relationship between two countries, from i to j , where w is 
the volume of the trade as provided by datasets fl 8 l [191. In 
what follow, we explain how the transition rules are assigned. 
In general, given a system of N nodes and a N x N matrix 
representing their interactions, the matrix is said balanced if 
if s° ut = s\ n , Vi £ {1,..., TV}, i.e., the sum of each col¬ 
umn is equal to the sum of its related row. Systems described 
by a balanced matrix can be considered isolated: nodes per¬ 
fectly balance their total in- and out-flows themselves, with¬ 
out needing any further exchange with the outside. However, 
most of the network representing real systems are not bal¬ 
anced. This means that s° ut s™ for at least one i such 
that i £ {1,..., N}. In these systems three different sets of 
nodes can be identified: the set of vertices with s° ut = s) n , 
the set of vertices with s° ut > s™ and the set of vertices with 

s out < s m 

In the first group, the total ingoing and outgoing 
fluxes of all nodes are balanced. In the second set, each vertex 
needs an extra incoming weight in order to balance in and out 
strengths. In this case, in order to balance the nodes in this set, 
we introduce an additional node/vertex, called source and la¬ 
beled with 0, providing extra-ingoing fluxes to all nodes in the 
set. The new link between the source and each node in the set 
will have the following weight: woi = s° ut — s" 1 . Similarly, 
a new vertex is introduced for the nodes in the third group in 
order to balance their in-strengths. It is called sink and labeled 
N + 1. The new links between this latter and the nodes in the 
set will have weights equal to: uujv+i) = s* n —s°“*. The ver¬ 
tices source and sink play the role of the internal economies 
of the countries. Trade surplus and trade deficit are absorbed 
by internal economies. For each node i with i £ {1,..., N} 
the total outgoing flow, v t , is given by: v, = s° ut + w^ N+ 1 ). 
Now, we introduce a N x N directed and weighted matrix M, 
such that rrijj = w t jVi. The elements of M, represent the one- 
step transition probability for a single particle (unit of value) 
to go from i to j. Higher powers of M express the transition 
probability from i to j in a given number of steps. We indicate 


with U^ s \ with S integer, the sum of the first S powers of M: 


s 

,(S)\ ^ 

< 2=0 


u (S) = « ; )i< i .,-<Ar = > 'Af« 


(/ -M s+1 ) 
(J — M) 


(A2) 


Since M is a sub-stochastic matrix, the series converges for 
S —» oc 1441 . The non-diagonal elements of U , ' S> represent 
the probabilities of reaching node j starting from node i within 
S stepstl0ll34|. Therefore, if we take into account all possi¬ 
ble paths of any length between two nodes (i.e. S = oc) 
eq.3 becomes = (J — M ) -1 . The element u[p en¬ 

ables us to compute the probability to come back to the source 
node i within S steps. We can therefore compute the cycling 
index of a vertex as the fraction of trade passing through a 
node/country i that returns statistically (directly or indirectly) 
to it within S steps, formally: 


(S) uu-1 iS ) 

11 ~ L t 


(A3) 


The global cycling index of the network (i.e. the fraction of 
traded value that returns, directly or indirectly, to a starting 
node) within S steps is given by: 


p(S') _ 


E, n 

J2 t v i 


E t if^ 

Ei«i 


(A4) 


(S) 

where 7 , = T^ v, represents the quantity of cyclic trade (mil¬ 
lions of current U.S. dollars), that returns to a node/country i 
within S steps. The value of T' fS! ranges between 0 and 1. 
The former case is observed when there is no trade that starts 
at some country i that comes back to it. The latter case repre¬ 
sent a systems where all trade come back to the starting coun¬ 
try. Note that the residual share of trade, up to 1, represent the 
share of trade that is acyclic, the starting country is different 
from the ending country. The quantity I ’' 00 - 1 has been used in 
ecological studies, with the aim of evaluating the total amount 
of cyclical matter in ecosystems lfl6l i23ll38l . 


4. Interpreting the cycling index. 

Ecosystems are open systems exchanging both matter and 
energy with a source and a sink. The above mentioned cy¬ 
cling index was developed to assess the share of matter that 
is recycled throughout food chains in an ecological network, 
from the primary producers (photosynthesis) to the top preda¬ 
tors and detritus feeders. Likewise, we can think of added- 
value as matter in food chains and look by means of cycling 
how this is conserved throughout the stages of production in¬ 
ternationally (GVCs), where the sink and the source of value 
are national economies. National economies play the role of 
sink and source in the parallel with ecosystems because the 
WTW, like ecological networks, is an open system: national 
economies play the role of sink and source in the parallel with 
ecosystems because the world trade network, like ecological 
networks, is an open system: matter is not conserved through 
every stage of the international production chain, whereby ev¬ 
ery country is a source of row material and a sink of waste. In 
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this paper the cycling index will be used to assess the GVCs. It 
is worth noting that the cycling index only measures, statisti¬ 
cally, the amount of value that returns to a starting point (coun¬ 
try) after a given number of steps (trading relationships), that 
is, it assesses the cyclical GVCs, from the markets of the raw 
material, to those of the final products. Three main method¬ 
ological approaches have been used to capture GVCs in the 
scientific literature: 1) international trade statistics on parts 
and components; 2) customs statistics on processing trade and 
3) international trade data combined with input-output (TO) 
tables. Amador and Cabral provides a detailed review of the 
existing methodologies @. More recently, network theory has 
been applied to disaggregate trade and 1-0 matrices to investi¬ 
gate GVCsl39ll43l. Markov chain theory has been previously 
applied to disaggregate trade to investigate the allometric scal¬ 
ing of networks and the structure of GVCsll4l l50l . Compared 
to the methods aimed at directly assessing the GVCs by mea¬ 
suring the traded value added in 1-0 matrices, our approach is 
different for the following reasons: 1) we do not measure just 
the share of incorporated value in exports/imports between 
pairs of countries, i.e. in bilateral relationships, we measure 
the value that is conserved throughout all possible paths in 
a network. 2) this is a statistical measure, thus, it does not 
rely on direct measures of value added (which can only be 
drawn by I/O matrices of countries). It has indeed the flaw 
of being a statistical measure of value, but it has the strength 
of assessing value through longer-than-one steps of trade (on 
the contrary, direct measure of value added can only assess 
one-step path, i.e. between two countries, inasmuch as any 
following step of the value added this becomes input of pro¬ 
duction). 3) it is a global (statistical) measure and asses all the 
possible paths of a given length. In other words, even when 
assessing two-step cycling, it is statistically relevant all the re¬ 
lationships of i and j with all other countries rather than just 
the relationship between country i and country j. 4) we sta¬ 
tistically asses (with cycling) the value conserved throughout 
a cycle, for all products, rather than one. This is a statistical 
measure based on aggregate trade because at every step of the 
production chain (generally, but not always) products change 
trade category, i.e. classification: iron, engines, cars, etc. To 
use an example from ecological networks - in a food chain, 
when we want to assess the amount of mass that is conveyed 
through one species (prey) to the other (predator) at every step 
of the chain, from primary producers (grass) to the last preda¬ 
tors (and decomposers), we cannot tag every atom and check 
every passage they make. We can only weight body mass of 
organisms through the food chain. If we know that species 
A feed 50% on species B and 50% on species C, we know 
that the atoms of the species A have 0.5 probability of coming 
from B and 0.5 of coming from C. We can do this for all the 
species of the food chain and we project this into a continu¬ 
ous, steady food relation, i.e., if species C feeds on species E 
for 50%, species A, even if it does not feed on species E, has 
0.25 probabilities of having atoms from species E. Upon this, 
we can calculate the probabilities of an atom to go from one 
species to the other through all the possible direct and indirect 
paths. This is referred to as transition matrix, and in the tran¬ 
sition matrix, we can calculate the share of atoms that make a 


cycle, i.e., that start from species A and come back to species 
A along all the possible paths (not only with species B and 
C, direct feeding, but also along species E, indirect feeding). 
Now, suppose we are not talking about atoms, but value of a 
product. If, for example, Italy sells cars to USA, where the en¬ 
gines of the Italian cars are produced, the share of value of car 
relative to engine is cyclical with USA. Suppose now that the 
USA buys iron from China and that Italy sells cars to China. 
Even if Italy does not buy directly iron from China, the share 
of the value of iron in the engine of the car is cyclical. 

5. Null model as a benchmark 

A popular and appropriate, to our scope, null model is 
the directed weighted configuration model indeed it preserves 
the observed intrinsic heterogeneity of vertices: all vertices 
have the same in-strength and out-strength as in the real 
network |[46i . In other words, this model preserves the in- 
and out-strength sequences separately, and, furthermore it pre¬ 
serves the total weight of the original network. To evaluate 
if our analysis if our analysis is sound and consistent. We 
proceed as follows: first, for each year we built the expected 
network of trade using the randomization approaches of the 
maximum-likelihood method also called exponential random 
graph, second, we evaluate fT') on the expected networks. 
See supporting information for the expected trend of T^T 
Note that in doing so we take into account the global trade 
growth and for each country the distribution of import and ex¬ 
port with foreign countries from 1960 to 201 l l46l . 

Appendix B: Supporting Information 
1. Methods 

In probability theory a Markov chain is a stochastic pro¬ 
cess defined on a discrete state space satisfying the Markov 
property. It is a set of random variables, representing the evo¬ 
lution of a certain system, without memory: each actual state 
of the system just depends on the previous one. The changes 
of states are called transitions and the Markov chain can be 
described by a Markov matrix, M, whose elements rn, :) rep¬ 
resent the probabilities from a state i to a state j (transition 
probabilities or transition rules). Given the memoryless na¬ 
ture of a Markov chain it is not possible to predict the state 
of a system, in a given time t, but it is possible to predict its 
statistical properties. 

In this paper the transition probabilities are set according to 
the elements of the matrix W that describes the World Trade 
Web, i.e. the trade network. A link, Wij, is assigned to any 
import/export relationship between two countries, where w is 
the volume of the trade. In what follow, we explain how the 
transition rules are assigned. 

In general, given a system of N nodes and a N x N matrix 
representing their interactions, the matrix is said balanced if if 
s° ut = s- n , V i € {1,..., N}, i.e., the sum of each column is 
equal to the sum of its related row. 
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Systems described by a balanced matrix can be considered 
isolated: nodes perfectly balance their total in- and out- flows 
themselves, without needing any further exchange with the 
outside. However, most of the network representing real sys¬ 
tems are not balanced. This means that s° ut ^ s™ for at least 
one i such that i £ {1,..., TV}. In these systems three dif¬ 
ferent sets of nodes can be identified: the set of vertices with 
s° ut = the set of vertices with s° ut > s™ and the set of 
vertices with s° ut < s 

In the first group, the total ingoing and outgoing fluxes of 
all nodes are balanced. In the second set, each vertex needs an 
extra incoming weight in order to balance in and out strengths. 
In this case, in order to balance the nodes in this set, we in¬ 
troduce an additional node/vertex, called source and labelled 
with 0, providing extra-ingoing fluxes to all nodes in the set. 
The new link between the source and each node in the set will 
have the following weight: woi = s° ut — s] n . Similarly, a new 
vertex is introduced for the nodes in the third group in order to 
balance their in-strengths. It is called sink and labelled TV + 1. 
The new links between this latter and the nodes in the set will 
have weights equal to: ru^jv+i) = s\ n — s° ut . 

For each node i with i £ {1,... ,7V} the total outgoing 
flow, Vi, is given by: Vi = s° ut +w^n+i)- Now, we introduce 
now a N x N directed and weighted matrix M, such that 
nriij = Wij/vi . We indicate with U^ s \ with S integer, the 
sum of the first S powers of M: 


(7 (s) = (ujf) i<i,j< N = I + M 1 + M 2 + ■ ■ ■ + M s = 


9=0 


{I-M s+1 ) 

( I-M) ' 

(Bl) 


Since M is a sub-stochastic matrix, the series converges for 
S oc 1441 . 

The diagonal elements of U < ' S ' ) allow the evaluation of the 
cycling index, Tj ; for each country according to the follow¬ 
ing r = u " (s) 1 . It is r | S ' 1 is a fraction of v%, therefore we 

can evaluate the total value of the flow (directly or indirectly) 
returning to a node as: 


7i 


(S) 


= rfV 


(B2) 


Let us assume a balanced and connected network. In net¬ 
work theory, two nodes i and j are strongly connected if there 
exists a path from i to j and vice-versa. This definition of the 
relation is reflexive, symmetric and transitive, dividing nodes 
in disjoint sets of equivalence. Each class is defined SCC and 
among them all nodes are strongly connected to each other, 
but no bidirectional path exists between the classes. 

First, we will prove by contradiction that a balanced and 
connected network can have one and only one SCC, i.e., each 
node is reachable from any other node. Thus, let us assume 
that we have a network, represented by the matrix W, with 
two distinct SCCs (see picture). 

We suppose that p of the N nodes of the networks belong to 
SCC\, while the remaining 2 vertices belong to SCC 2 - The 
network is connected, thus exists at least a link w l f ; from a 
node i £ SCC± to a node k £ SCC 2 - 

Since that the network is balanced we can write the sum of 
each row of matrix W equal to the sum of the related column. 
We write down this equality for the nodes in SCC\, and we 
explicit the weights related to the node i: 


p 

y^w-u + w u = 

j =i 

p 

Y, W-2;, +W 2 i = 

3 =1 
V 

y w pj + w pi = 

3=1 


P 

y Wji+wn 

3=1 

p 

y Wj 2 +w i2 

3=1 

p 

y Wj p + w ip . 
3=1 


If we sum up all the left side members and right side mem¬ 
bers in the equalities and simplify, we obtain: Wji = 

Ej Wij- We observe that V'' ir ,, = s} n , while Ej w v = 
s° ut — Wik■ The network is balanced by hypothesis, therefore 
Wik = 0. This result contradicts the assumption that the net¬ 
work is connected, leading to the conclusion that it is impos¬ 
sible to have more than one SCC in a balanced and connected 
network. 


The global cycling index within S steps is given by T^ = 


a. Two extreme cases 

The quantity has a lower bound, 0, and upper bound, 
1. The two extreme cases are = 0 and lim = 1. 

S—too 

The former is observed in Directed Acyclic Graphs (DAG), 
and the latter in balanced networks. While the case = 0 
is straightforward due to the fact that circular paths do not 
exist in a DAG, the second case is not trivial to picture. We 
give a quite intuitive idea of the proof in two steps. 


The second step of the proof consists in showing that for 
a balanced network u u —> oo for S —> oo. We know that 
a balanced and connected network is strongly connected and 
implies that the related Markov chain is irreducible, in other 
terms it is a single communicating class. A state i of a Markov 
chain is called recurrent if the probability to start from i and 
come back to it is equal to 1. In particular, if the chain is finite 
and irreducible, then all states are recurrent [? ]. 


In an irreducible Markov chain with finite number of nodes 
the nodes are recurrent. N odes are recurrent if and only if 

Es 


^°° ^ = oo, hence eq. 


B2 


i _ ( becomes: T^ 00 ) = 

1 yy = 1- We conclude that for a balanced connected net¬ 
work r(°°) = l. 
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(b) 

FIG. 5. From 1960 to 2012, the trends of F^ 2 ' (orange line), T® 
(red line) and r^ 4) (purple line) for the USA (panel a) and Germany 
(panel b). 


Country 

p(2) 

L i 

p(3) 

L i 

W4) 

L i 

p(°°) 

1 i 

CHI (1961) 

0.03(0.8) 

0.08(0.5) 

0.10(0.4) 

0.18(0.2) 

FRA (1961) 

0.63 

0.53 

0.47 

0.18(0.2) 

GER (1971) 

0.25(0.1) 

0.26(0.1) 

0.26(0.1) 

0.26(0.1) 

JAP (1961) 

0.76 

0.75 

0.75 

0.44 

U.K. (1961) 

0.16(0.2) 

0.19(0.2) 

0.2(0.1) 

0.22(0.1) 

RUS (1990) 

-0.57(0.02) 

-0.54 (0.01) 

-0.53(0.01) 

-0.54(0.01) 

USA (1961) 

0.31(0.02) 

0.31(0.03) 

0.31(0.03) 

0.28(0.04) 


TABLE II. For different countries, the table reports the Pearson’s 
correlation index between values of F - r) (for r = 2, 3, 4, oo) and 
GDP growth at constant price, in parenthesis the p-value (95% con¬ 
fidence, and p-value lower than or equal to 10' 5 are not reported) . 
GDP growth data are from World Bank; the starting year of World 
Bank statistic is reported near each country. 


2. Supplementary Discussion 


The method described to build M has been used in eco¬ 
logical studies, with the aim of evaluating the total amount of 
matter cycling in the ecosystem with the quantity r (oc) [? ? 
? ]. Recently this approach has been employed to detect the 
presence of allometry patterns both in ecosystems and World 


Trade Web (WTW)[? ? ]. The novelty of our methodology 
stands in considering all possible orders of cycling paths, with 
S to {2, oo}, in order to evaluate the Global Chain Value, 
and the share of cyclic trade in the world economy. 

In this perspective, the three different sets of nodes defined 
can be thought as groups of countries with a perfect balance 
between export and import ( s° ut = s'" 1 ), a surplus of exports 
s° ut > sl n or of imports s° ut < s" 1 . According to this inter¬ 
pretation, the sink and source nodes just represent the total of 
the World’s national economies. 

In the paper we showed the strong correlation between the 
global cycling index and the crude oil price. This is the main 
result of our work in terms of policy implications, nevertheless 
the single country analysis reveals interesting information. 

The two quantities 7 ] ’ and T ■ ' can be used to evaluate 
the performances of each country from 1960 to 2011. Their 
trends can reveal significant differences in the structure of the 
chain values and shed light onto the distances in the underly¬ 
ing economic system. 

(S) 

Figure 4 shows the fluctuations of a sample of T) for USA 
(left panel) and Germany (right panel). 

For both countries, a big share of exports returns directly 
or indirectly to the country, however the plots depict different 
cycling patterns. Specifically, taking into account ciclyc paths 
longer than 2 steps, USA seems characterized by a longer 
chain of production then Germany. Indeed, the share of USA 
cyclic paths equal to 4 steps is bigger than the 3 steps one, 
while the majority of German cyclic paths show length equal 
to 3. Furthermore the trends of 1’,; curves show a different re¬ 
sponse of the two countries to the 2008 crisis. A decrease of 
the share of cyclic flows is observable since the end of 90s and 
it culminates into a dramatic drop in 2008. 

The overall percentage variation is around —15%. On the 
other hand, the German declining trend starts around 1990 and 
seems to stabilize from 2000 on, and the 2008 crisis appears 
negligible. The overall percentage variation of the share of 
cycling paths is around —5%. 

We could argue that the cycling coefficients and the growth 
of Gross Domestic Product of a country are related. How¬ 
ever this connection strongly varies from country to country 
without a clear pattern. Indeed, in table (II) we show the 

Pearson’s correlation coefficient between the values of GDP 

(S) 

growth (p-values in brackets) at constant price and U , with 
S £ {2, 3,4, 00 }, for six countries: China, France, Germany, 
Japan, United Kingdom, Russia and USA (GDP data comes 
from World Bank, www.worldbank.org). Most of the coun¬ 
tries do not show any significant correlation or just a slight 
one; France and Japan exhibit a strong positive correlation be¬ 
tween the two variables, while for Russia we found a quite 
high negative correlation. It means that the information cap¬ 
tured by gamma substantially differs from what can be drawn 

from the analysis of the GDP growth. 

(S) (S) 

The 7 } and I; can give information about the perfor¬ 
mance of single countries from 1960 to 2011. It is possible 
to rank countries according to their values. In figure (5) we 
show this rank, for S £ { 2 , inf}, of the same seven coun¬ 
tries: China, Germany, France, UK, Japan, Russia, and USA. 

We chose these countries for their interesting patterns in 
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(0 


(d) 


to the value of 7 ® (panel a), (panel b), F^ 1 2 ' (panel c) and r^ 00 ^ (panel d). 


FIG. 6 . The ranking position for seven countries according 
The data range covers from 1960 to 2011. 

• (S) (S) 

order to reveal how the variables 7 ? and IT ; can explain 
structural changes in countries. Russia position is always 
higher in the relative than in the absolute rank, and for S = 2 
than for S = 00. Moreover we can observe a decline in both 
starting in late 80s and a rise since around 2000. Results seem 
in line with the USSR dissolution, explained by a Russian 
chain value made up substantially by bilateral exchanges with 
the other Soviet Republics. 

China is characterized by a remarkable increasing trend in 


the last decades, reaching the top of the rank for ''/ i ° c ' > in the 
last years of our analysis. In the same period, Germany over¬ 
takes USA share of cyclic trade (S = 00), as emerges from 
fig. 5 panel d. 

We remark on the performance of China and Germany. The 
former has been on the rise in the last decades and we point out 
the pervasiveness of its chain value in absolute value ( 7 ^°°^, 
fig. 5 panel b). In the last year regarding S = 00, Germany 
has outdone USA for share of cyclic trade (fig. 5 panel d). 
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